Random sequential adsorption of shrinking or spreading particles 
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We present a model of one-dimensional irreversible adsorption in which particles once adsorbed 
immediately shrink to a smaller size or expand to a larger size. Exact solutions for the fill factor and 
the particle number variance as a function of the size change are obtained. Results are compared 
with approximate analytical solutions. 
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I. INTRODUCTION 

Random sequential adsorption (RSA) is an attractive 
model for a number of physical phenomena, including 
such different application as information processing [1| 
and particle branching in impact ionization [2j ■ The sim- 
plest example of RSA is the car parking problem (CPP). 
Of interest is the average number of particles ("cars") 
adsorbed on a long line, as well as the variance of this 
number (see Refs. li for the review). In what follows, 
we shall use the term standard RSA for the classical 1- 
dimensional model corresponding to particles of fixed size 
that arrive randomly on a line and are deposited if empty 
space is available and rejected otherwise. Extensions of 
this model include RSA with particles expanding in the 
adsorption process [U , two-size particle adsorption [f| 0] , 
and also RSA with an arbitrary particle-size distribution 
function [8[. 

Few of these models, especially physically relevant 
ones, have an exact solution. As a rule, only the fill 
factor is determined. However, for a number of appli- 
cations fluctuations are of major importance. An exam- 
ple of such an application is the very important practi- 
cal problem of particle energy branching (PEB) where 
high-energy particle propagates in an absorbing medium 
and multiplies producing secondary electron-hole (c-h) 
pairs. The energy distribution of secondary particles is 
random to a good approximation. Multiplication pro- 
ceeds so long as the particle energy is above the impact 
ionization threshold [9(. This connection was noticed as 
early as in 1965 by W. van Roosbroek The exten- 
sion of RSA model proposed in , known as the " crazy 
carpenter model" was further exploited in 0, E| • 

The PEB process can be considered in terms of a CPP 
if one identifies the initial particle kinetic energy with an 
available parking length and the pair creation energy as 
the car size. Full equivalence to CPP requires further 
that only one of the secondary particles takes on signifi- 
cant energy, otherwise one has to consider simultaneous 
random parking of two cars in one event [13 ]. 

The number of created electron-hole pairs N in PEB 
serves to evaluate the initial energy. Variance of this 
number limits the accuracy of energy measurements. 



Both the yield N and the c-h pair variance var(iV) = 

(N — N) 2 are proportional to initial energy. The ratio of 
the e-h pair variance to the yield is called the Fano fac- 
tor [l3[ . The Fano factor $ (for a Poisson-distributed N 
one has $ = 1) is a parameter that quantifies the energy 
resolution of high-energy particle detectors. For semicon- 
ductor crystals, the PEB problem has many additional 
complications due to phonon losses, as well as features in 
energy dependencies of the particle density of states and 
the impact ionization matrix element. 

Earlier attempts to evaluate the Fano factor for the 
PEB problem in semiconductors emplo yed widely differ- 
ent approaches (compare @, [n| and [111)- To obtain 
agreement with the experimentally observed $ w 0.1 (for 
semiconductors), different fairly rude and unjustified as- 
sumptions were made, so that the numerical coincidence 
is of little value. Also available are numerical calculations 
pdl . [TBI . However, the relative importance of various fac- 
tors (e.g. phonon contribution to $) remains question- 
able within the numerical models, while the precision of 
results is difficult to assess. Evaluation of the Fano fac- 
tor is important for predicting the energy resolution of 
detectors, especially those based on new materials and 
new principles [l6j| . 

An important aspect of the impact ionization process 
is the fact that the threshold impact ionization energy is 
usually larger than the minimum energy E g needed for 
c-h pair creation. This difference arises from kinematic 
restrictions imposed by momentum conservation, so that, 
e.g., for equal effective masses of electrons and holes the 
minimal energy of pair production is 3/2E g . In crystals 
with non-equal electron and hole masses the threshold 
can vary from E g to 2E g . In terms of the RSA problem 
this is equivalent to the particle shrinking immediately 
after the adsorption (parking), making larger available 
length for the next adsorption event. 

For the PEB problem only particle shrinking is per- 
tinent. However, for a general RSA problem both the 
shrinking and the expansion of particles are relevant, due 
to such factors as, e.g., repulsion between particles and 
surface attraction. 

In this paper, we use a recursive approach to consider 
the general RSA problem for particles that either shrink 
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or expand immediately upon adsorption. We present an 
analytical solution for both the average filling factor and 
its variance. The results are compared with approximate 
solutions obtained by the methods employed in earlier 
approaches to the energy branching problem 0, 

MM 

Il4| and, for expanding particles, also with the results 
of a kinetic approach used in Ref. Exact solutions 
presented here allow to assess the errors introduced by 
the adopted approximations in the evaluation of both the 
yield and the variance. 



II. STATISTICS OF ADSORPTION 

In this section we develop an analytical model of ad- 
sorption for shrinking or expanding particles. We shall 
be using the recursive technique [3]. Besides being the 
most direct approach, this technique has the further ad- 
vantage of being applicable for the evaluation of both 
the average filling factor and the variance. The recur- 
sive approach has been confirmed by direct Monte-Carlo 
computer simulations liol. [Ti| for similar problems. 

Consider adsorption of a particle on an initially empty 
line of progressively growing length x. The final (after 
shrinking or expansion) size of the adsorbed particle is 
taken equal unity. Denote by w x the random variable 
corresponding to the total wasted length for some config- 
uration of adsorbed particles in the jamming limit, when 
all gaps are smaller than the minimal length needed for a 
particle to be adsorbed. In the process of sequential ad- 
sorption deposition of additional particle generates new 
gaps with the same distribution of the gap size. Equa- 
tions for the expected values of the moments of waste dis- 
tribution can be obtained usin g mo ment generation func- 
tion, ^(A, w x ) = Eexp(Aw 2: ) [l7[. Moments of w x are 
obtained by differentiating with respect to A. When 
a particle is adsorbed into gap x, two new gaps appear, 
the size of the gaps being y and x — 1 — y. Therefore for 

we have 



(1) 



In this equation, the angular brackets denote averaging 
over the distribution of y, characterized by a one-particle 
gap distribution function (OPDF), p(y\x), that equals 
to the probability density of creating an interval y in a 
one-particle deposition into the line x, viz. 



< $(\,w y ) >= 



y{\,w y )p(y\x)dy. (2) 



The randomness in y is induced by a two-step random 
process, (i) the random choice of the adsorption coordi- 
nate and (ii) the subsequent size change, which depends 
on the adopted model of shrinkage/expansion. 

Using Eq. (JTJ) in calculating the first derivative of 
ty(\,w x ) at A = 0, we find the first moment (the mean 
value) of waste, ~EiW x = w(x), in the form 

w(x) =< w(y) > + < w(x - 1 - y) >= 2 < w(y) >, (3) 



y =0-1) 

J mm, r y f 



(a) 







y = x-l 

J max, r 


1 










/ 



p(y\x) 




FIG. 1: Illustration (a) of particle adsorption near the edge of 
a gap of size x in the model where the particle is assumed to 
shrink upon adsorption (from size I to unit size) by randomly 
retracting one of its ends (the retracted part is marked by the 
dark shading, the final state is shown as a white brick of unit 
length); y m in,r and y m ax,r are minimal and maximal gaps 
created by particles adsorbed at left end of the length x and 
retracting to right endpoint ; part (b) shows the one particle 
gap probability distribution function or OPDF, p(y\x). 



the latter equation being due to symmetry. The second 
moment, u(x) = Ew^, needed for evaluation of variance, 
is obtained by taking the second derivative, 



u{x) = 2 < u(y) > +2 < w(y)w(x — 1 — y) > 



(4) 



Equivalent equations can be derived for the mean filled 
length f(x), which for a unit particle size equals the num- 
ber of particles adsorbed on the line of length x. The ran- 
dom variables w x and f x are complementary in the sense 
that w x = x — f x and, therefore, equations for f(x) can 
be obtained from Eqs. and Q by the substitution 
w(x) = x — f(x). 

For the standard RSA problem the OPDF is a homo- 
geneous distribution of y in the interval {0, x — 1| ). The 
corresponding equations were obtained in [l|, [lof in the 
context of RSA and in 0, H3] for the energy branching 
problem. 

An alternative approach to RSA is to consider a ki- 
netic (or rate) equation that describes the sequential de- 
position of particles. In the kinetic approach the gap 
size distribution function G(x,t) representing the aver- 
age density of voids of the length between x and x + dx 
at a time t obeys the equation fl8j 



dG{x,t) 



dt 



G(x,t) I dsp(s) / dyF{y,x-y - s\s) 
o Jo 

)f dsp(s)F(x,y-x-s\s). (5) 



Here p(s) is the adsorbed particle distribution function, 
and F(y, x — y — s\s) is the deposition probability, that 



determines the average rate, 



R(y\x) 



dsp(s)F(y,x - y - s\s), 



(G) 



at which the initial length x is destroyed by the deposi- 
tion of a particle producing a gap y. One can see that for 
a fixed x = xo the rate R and our OPDF are proportional 
to each other, p(y\xo) oc R(y\xo), the only difference be- 
ing due to the fact that p(y\xo) is normalized to unity. 

For the case when the randomly adsorbed particles 
change their size, the OPDF depends on the particu- 
lar model of size transformation in the adsorbed state. 
Several such models are discussed below. 

Shrinking particles. First, we consider a model in 
which particles of initial length I > 1 shrink to the length 
of unity by randomly retracting one of their endpoints 
(cither on the left or on the right with equal probability) . 

In this model, no particles are adsorbed for x < I, i.e., 
small intervals are wasted entirely, w{x) = x. When the 
interval length x reaches I, then two gaps of size (I — l)/2 
are created. For an arbitrary x > I, adsorption of a 
particle creates two new gaps y and x — 1 — y. Note 
that the probability of creating gaps near the edges of the 
initial interval (i.e., for < y < I — 1 and x — I < y < x— 1) 
is twice smaller, since they get contributions only from 
particles that retract the right endpoint (at the right side 
of the interval a;) and the left endpoint (at the left side 
of x), respectively, see Fig. 1. This model can be used 
for the PEB problem, to account for the effect of lower 
density of states at low particle energies. 

For sufficiently large intervals, x > 21 — 1, the mean 
waste is described by the following equation 



w{x) 



I 



- I 



x — l 



Z-l 



w(y)dy 



' X — l 



(7) 



Note that Eq. ([7]) does not work in the region I < x < 
21 — I, where the gap distribution p{y\x) is simply uniform 
within y < x — I and I — 1 < y < x — 1. However, since 
the adsorption of the second particle starts at x > I + 1, 
knowledge of the variation of w{x) at small x < I + 1 
is sufficient to obtain an exact solution of Eq. ([7]) — so 
long as 21 — 1 < I + 1 (or I < 2). In what follows we shall 
confine ourselves to the case Z < 2. 

According to (0| the equation for the second moment 
u(x) is obtained from Eq. ([7]) by the substitution w{y) — > 
u{y) + w(y)w(x - l-y). 

Another possible shrinking model results when the ad- 
sorbed particle is assumed to shrink symmetrically about 
its center. In this case the OPDF remains constant 
within the interval (I - l)/2 < y < x - (I + l)/2. This 
model accounts for the particle repulsion in a model of 
hard disks, when the distance between the particle cen- 
ters can not be smaller than a certain length that serves 
as their effective size. That length can often be larger 
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FIG. 2: Adsorption of expanding particles with the transfor- 
mation rules corresponding to a symmetric expansion a =>• 1 
when size permits and near an edge an asymmetric expansion 
to fill the available space; (a) illustration of particle adsorp- 
tion near the edge of a gap of size x; (b) one-particle gap 
distribution function, p(y\x), shows 5-function singularities 
at y = x — 1 and y = 0. 



than the actual particle size used in the description of 
the resulting coverage. A similar model reasonably de- 
scribes the PEB process in semiconductors, when both 
secondary particles produced by impact ionization are of 
the same mass and acquire equal kinetic energies. 
With this OPDF, equation ^ takes the form 



w(x) 



x — I 



x-{l+l)/2 



w(y)dy. 



(8) 



'(l-l)/2 

Similarly, equation for the second moment u(x) reads 

r x-(l+l)/2 



u(x) 



1 



u(y)dy 

(l-l)/2 

2 r x-(i+i)/2 



x — I 



(l-l)/2 



w{y)w(x - I- y)dy. (9) 



These equations are different from those obtained for 
the standard RSA problem. The modification arises due 
to the edge effect and is important at the jamming limit, 
when gaps are minimal. For particles shrinking to the 
center, the edge effect results in stronger restrictions on 
the gap size. 

Expanding particles. Here we consider a model in 
which particles of initial length a < I given sufficient 
space expand symmetrically to unit length. If the particle 
is placed in a gap of size a < x < 1, it fills it completely 
If it is placed in a gap x > 1 near its edge, it expands 
asymmetrically to unit length. This model represents size 
transformation due to an attractive force from the surface 
for the case when interaction between the particles is 
negligible. 

To define the ODPF for this model, we note that given 
the initial length x, coordinates of the centers of ad- 
sorbed particles are homogeneously distributed within 
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a/2, x — a/2. For particles adsorbed not too close to the 
gap edges, the gap distribution function remains uniform 
in the interval < y < x — 1. However, all adsorbed 
particles whose centers fall within an interval (a/2,1/2) 
from an edge will produce a gap of the same size x — 1, 
see Fig. 2, as well as a "gap" of zero width. This results 
in a singularity at y = x — 1 (and another one at y = 0) 
in the gap distribution function (existence of this singu- 
larity for expanding particles was already noted in [5|). 
With this ODPF, Eq. (J3j) acquires the form 



w{x) 



x — a 



w(y)dy 



1 



-w(x-l). (10) 



Similarly, the second moment equation becomes 



u(x) 



x — a j 
2 



x — a 



1 — a 

uyyjdy H u(x — 1) 

x — a 



w(y)w(x - 1 - y)dy. (11) 



Note that the particle expansion and edge effect do not 
influence the second term of Eq. (fTTj), since the product 
w(y)w(x — 1 — y) = vanishes near the edges. Neverthe- 
less, we note that both Eqs. (jTD|l and (fTTj) are modified 
by the edge effect. 

In all models of size transformation considered above 
the ODPF remains constant within certain intervals. It 
is this feature that enables exact solution of the recursion 
equations. 

It is worthwhile to stress that due to the self-averaging 
nature of the filled length (and the waste length) in the 
limit x — > oo solution of the averaged (hence approxi- 
mate) recursion equations gives exact results. The same 
is true for the kinetic approach in the limit t — ► oo. Both 
approaches are, therefore, equivalent for the calculation 
of the fill factor. 



III. FILLING FACTOR 

We evaluate the filling factor in the so-called "jamming 
limit" — corresponding to the situation when every gap 
capable of adsorbing a particle has done so. The average 
density of particles saturates in the jamming limit. 

Shrinking particles. We shall first focus on the 
model of particles shrinking by randomly retracting one 
of the endpoints. The case of particles shrinking to the 
center will be discussed subsequently. We have to solve 
Eq. with appropriate boundary conditions that orig- 
inate from the region of small available adsorption inter- 
val x. Consider an initially empty line of progressively 
growing length x. For x < I the entire interval is wasted, 
w(x) = x. For I < x < I + 1 only one particle will be 
adsorbed producing waste equal to w(x) = x— 1. Clearly, 
for x < I + 1 the waste is fixed and does not fluctuate. 
For x > I + 1 the probability that two particles will be 
sequentially adsorbed grows steadily and so does the av- 
erage waste, which also begins to fluctuate. At large x 




FIG. 3: Average filled length / and the wasted length w for 
shrinking particles as functions of the length x of the adsorp- 
tion interval assumed initially empty. The assumed shrinking 
model is particles randomly retracting their endpoints. The 
results are obtained by iterating Eq. (9) with the assumed 
shrinking ratio of I — 1.6. 



both the wasted length and the average covered length 
grow linearly with x, so that the average waste per par- 
ticle or per unit adsorption interval remains constant. 

For particles whose final size is fixed it is natural to 
study the mean number of particles f(x) adsorbed in a 
line of length x. A convenient equation for f(x) valid for 
x > 1 is obtained by substituting w = x — f(x) in Eq. 
([7]) and making the replacement x — > x + I: 



f(x + l) = l + 



2 

x Ji 



c{x.l} 



1 



x+l—l 



X J rnax{x,l} 



f(y) + - I f(y)dy- 

(12) 

In deriving Eq. (fT2]) use has been made of the initial 
condition 



fO, 0<x<l 
I{) IT l<x<l + l 



(13) 



Evaluation of f(x) at small x > 1 + I is readily done by 
a repeated iteration procedure, going from the small to 
progressively larger length sizes. Results of the numerical 
recursion are shown in Fig. 3. Note that already at x rs 7 
the variations of /(x) and w(x) are very close to linear. 

An exact solution of Eq. (fl~2|) , which will be used to 
calculate variances, can be obtained using Laplace trans- 
formation. Multiplying Eq. l[T2"|) by x, taking the Laplace 
transform and using the initial conditions specified by Eq. 
(|13[) . we obtain an equation of the form 

- 4- {e pl F{p)) = ~ 2 + ~F(p) + -e^F(p) (14) 
dp p z p p 

for the Laplace transform of /(x), 

/•OO 

F(p) = / e-v x f{x)dx. (15) 
Jo 
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Rearranging terms and multiplying by e pl , Eq. (fT4)) can 
be rewritten in the form 



F'{p) 



I + - (e- pl + e- p ) 
P 



exp(-pZ) 

F yp) = -5 — • ( 16 ) 



The solution of Eq. I|16[) satisfying the boundary condi- 
tion at p — > co, 



^(P)| P - 



1 

— e 



-p i 



(17) 



which follows from the known variation of f{x) at small 
x, can be obtained in a straightforward manner: 



F(p) 



exp(-p I) 

' p 2 Mp) 



/3 r (u)du, 



(18) 



where 



P r {p) = exp 



2 — exp(— w) — exp(-vl) 



dv 



(19) 

In the analysis of the solution we follow the approach 
of Ref. [H, based on the use of Karamata's Tauberian 
theorem, see e.g. [lj|, p. 37, for the asymptotic growth 
rate of steadily growing functions. According to that the- 
orem, in order to obtain the asymptotic behavior of the 
filled length, f(x) (or the gap, w(x), or of the variances of 
these functions) it is sufficient to have a Laurent power- 
series expansion at small p of the Laplace transforms of 
these functions (possibly cut at small x by a Heaviside 
step- function factor). Further mathematical details of 
this type of analysis can be found in [l[ . 

Function F(p) is analytic at all p ^ and has a second- 
order pole at p — with the following expansion as p — > 0: 



F(p) 



where 



I3 r (p)dp. 



(21) 



To calculate f(x) at large x, we take the inverse Laplace 
transformation of the asymptotic expansion (|20p . This 
gives 

f(x) = a ffi x + otf >0 - 1, (22) 
with an exponentially small error term. Whence we have 



w(x) 



•),o( x + 1)> a w,a = 1 - a/,o- 



(23) 



In the limit 1 = 1, equation (|2"Tj) gives the so-called jam- 
ming filling factor R for the standard RSA, otffi{l = 1) = 
R = 0.74759 •• • (also called the Renyi constant [20l|). 
The filling factor (saturation coverage), calculated with 
Eq. (|2ip as a function of the shrinkage ratio I, is shown 
in Fig. 4 (curve 2). As expected, the shrinking causes 
a decrease of the filling factor with /. Owing to the fact 



average gap approximation 




1.2 1.4 1.6 1.: 

Shrinking ratio 

FIG. 4: Filling factor a/,o(i) as function of the shrinking 
ratio / for different shrinking models. Curve (1) describes 
particles shrinking to their centers. For the model of particles 
randomly retracting their endpoints curve (2) shows the exact 
solution and curve (3) the linear approximation. The upper 
curve corresponds to the average-gap approximation. 



that the correction terms are exponentially small, our 
asymptotic solutions (|22|) and (f23|) are extremely close 
(within less than 2 • 10~ 4 for x > 7) to the exact solution 
obtained by direct recursion (Fig. 3). This accounts for 
the linearity of }{x) and w{x) at high x evident in Fig. 
3. 

For particles shrinking to their centers, we use Eq. ([8|) 
and obtain 



/(*) 



- I 



:-l 



f[y + 



i - i 



dy. 



(24) 



(20) Equation (f2~4"|) can also be written in the form 



f x + l + 



l + l 



x+(l + l)/2 h 



x+l 



f (y) dy- 



(25) 

Equation (|2"5|) exhibits a recursion period (Z + l)/2 which 
is to be compared to a recursion period of unity for the 
standard RSA. The initial condition for Eq. (|2~5"]) is given 
by Eq. (|13|) and the resultant solution is similar in form 
to Eq. (fTS)) with the following replacement: f3 r (p) —> 
/3 c (p), where 



Pc (p) = cxp 



From Eqs. 





and (|21[) it follows that 
2 



a f,o 



l + l 



R. 



(26) 



(27) 



Equation (|2T)l represents the exact solution for this model 
in a simple analytical form. The fill factor calculated 
with Eq. (|2T|) is shown in Fig. 4 as curve 1. As can 
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be expected, the shrinking effect is smaller for particles 
shrinking by randomly retracting their endpoints than 
for particles shrinking to their centers. This is due to the 
reduced contribution of edge regions. 

As noted above, the model of particle shrinking to the 
center accounts for a strong repulsion between the parti- 
cles acting as hard disks. The minimal distance between 
the particle centers (an effective particle size) equals to 
(I + l)/2. Therefore the number of a unit size particles 
adsorbed in a line with fixed length and the resulting 
jamming limit fill factor are reduced by 2/(1 + 1). 

An alternative approximate way to obtain an asymp- 
totic solution of Eqs. (fT2"| . (|25[) is to seek the solution in 
the linear form f(x) = ax + b in the region x > nl + 1, 
for some n, while using a function obtained by direct it- 
erations for x < nl + 1. This recipe also gives a very 
good result starting from n = 1. Thus, say, for parti- 
cles shrinking to centers we obtain af^ = 2/(1 + l)R ( - 1 \ 

with i?W = 0.75, quite close to the exact Renyi constant. 
A less accurate approach is to merely match the linear 
asymptotic to the recursion result for x < nl + 1. In this 
case, one would need n — 2 to obtain a similar precision. 

In calculations of the quantum yield and the Fano fac- 
tor for the PEB problem, the complexity of equations 
often inspires even more radical approximations, based 
on estimates of the average losses, see e.g. @. In terms 
of the RSA problem, this is equivalent to assuming the 
filling factor a^o in the form afo = R(R + a w ,o) with 
the wasted length growing linearly with Z, i.e. taking 
ct w fi = (I — R) ■ I. Comparing the exact and the approx- 
imate solutions of Eq. (f2"4")) (see Fig. 4), we see that this 
approach severely underestimates the effect of shrinking. 

Expanding particles. We consider now the case 
when the initial particle size is a < 1. If the interval 
at the adsorbing line is a < x < 1, the particle fills 
it entirely (obviously, without fluctuations). If x > 1 
the particle expands to the length of unity. Evidently, 
for x < a the whole length is wasted, w(x) = x. For 
a < x < 1 no waste occurs, while for 1 < x < 1 + a 
the waste equals w = x — 1 and does not fluctuate. For 
x > a + 1 the probability of two particles being adsorbed 
in this space steadily grows with x and so does the aver- 
age waste. Thus, the initial conditions to Eq. (fit)]) are of 
the form 



w(x) 




< x < a 
a < x < 1 
1 < x < 1 + a 



(28) 



Fig. 5 shows the functions w(x) and f(x), obtained 
by direct iterations of Eq. (flU)) . Note a step-like feature 
in f(x) and w(x) at x = 1 + a (which is replicated at 
x = n + a with ever smaller amplitude) and a decrease 
of the w(x) in the region x > 1 + a. These are "recursive 
replicas" of the singularity in OPDF discussed above and 
the gap in w(x) at a < x < 1 . 

To calculate the variance we shall need an exact so- 
lution for w(x). It can be obtained by taking Laplace 




FIG. 5: Average filled length / and the wasted length w for 
expanding particles as functions of the length x of the adsorp- 
tion interval assumed initially empty. Results are evaluated 
by recursion with Eq. (8) for an assumed expansion ratio of 
I = 1/0.6. 



transformation of Eq. (fT0|) . First, it is convenient to 
multiply Eq. (|11| by x — a and make the substitution 
x — > x + 1. Taking Laplace transformation, we find that 
the function 



/oo 
e- px w{x)dx, 



(29) 



(which is the Laplace transform of w(x) cut at small 
x < 1 by multiplying with a step function) satisfies the 
following equation 



— - +l-a 
dp 



e pl F w (p)) = -F w {p) 



+(l-a)F w (p)+[- + l-a)J{jp), (30) 



where 



J(P) 



w(x)e 



-px 



(31) 



Using Eq. (|28|) to calculate J{p), we can re- write Eq. 
(|3"P|) in the form 



K(p) 

where 
and 



a + —b{ap)\ F W (P) = - C M^-G w (p), 



G w {p) = -b(ap)J 1 (ap), 
P 



(32) 
(33) 



1 f ap 

b(ap) = 1 + - (1 - a)p, Ji(ap) = / te^dt. (34) 

2 Jo 
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for the RSA problem with restricted particle expansion 
(their model allows expansion only when the required gap 
is fully available, see Appendix). As can be expected, 
the restricted expansion gives a smaller increase of the 
fill factor. 

Figure 6 also shows (curve 3) the approximate results 
obtained by neglecting edge effects; in this approxima- 
tion the expansion effect is reduced. Much better results 
are obtained in the approximation (shown by curve 2) 
obtained by using as a solution of Eq. (TTU)) a linear de- 
pendence at x > n + a recursively combined with the 
exact solution for x < n + a. Good agreement with the 
exact result is obtained already for n = 1. 



1.4 1.6 1.8 

Particle size ratio 



FIG. 6: Filling factor q/.o as a function of the expansion 
ratio / = 1/a; curve (1) shows the exact result evaluated with 
Eq. (34), curve (2) the approximate solution of Eq. (8), 
obtained by linear-to-iterative matching, curve (3) describes 
the approximate solution neglecting edge effects. Curve (4) 
corresponds to the restricted expansion model of Ref. [|| (see 
also Appendix). 



Solution of Eq 

f w (p) = 

where 

P(p) = exp 



is of the form 
exp(-p) 



p 2 (3(p) 



(3{u)G w {u)du, (35) 



b(av)dv 



(36) 



Function F w (p) is analytic at all p ^ and has a second- 
order pole at p = 0. As p — > 0, the following asymptotic 
expansion holds: 



where now 



F w {p) = 



O-wfi ot w ,o 



p 



0( P ), 



(3{p)G w {p)dp. 



Applying the inverse Laplace transformation to Eq. (|3~7| 
we bring w(x) at large x into the form w(x) = a W fl(x + l) 
with an exponentially small error term. The filling factor 
a/.o = 1 — a w fi in the jamming limit is then given by 



a f,o 



hp)b(ap) 



-ap 



+ ae- ap dp (39) 



P 



again reduces to the Renyi con- 



For a = 1, equation 
stant R. 

The filling factor a f y o calculated with Eq. (f39|) is plot- 
ted in Fig. 6 as function of the expansion ratio / = 1/a. 
We see that expansion causes an increase of the filling fac- 
tor with I, Also shown are the results of Boycr et al. Q 



IV. VARIANCE AND FANO FACTOR 

Shrinking particles. We shall consider in detail the 
model of particles shrinking by randomly retracting their 
endpoints, and then briefly discuss the case of shrinking 
to the center. 

Equation ^ can be readily transformed into an equa- 
tion for the expected value of the occupied length squared 
v = E/ 2 . Master equation for v(x) is of the form 



v(x) 



x — I 



x — I 



[v(y) + 2f(x)]dy 



x-l 



x-l 



[v(y) + 2f(x)]dy 



x — I 



x — l 



(40) 



f(y)f{y -x- l)dy. 



In deriving Eq. (j40|) we used the initial conditions ([13 
The Laplace transform M(p) = L(v(x)) satisfies 



(37) Af'(p) 



l + -(e 
P 



- P i 



M(p) 



exp(-pZ) 



where 



Rfip), 

(41) 



(38) Rf(p) = l+2 P F(p) (e^- 1 ) + l)+2p 2 F 2 {p)e 



p(i-i) 



(42) 



with F(p) defined by Eq. (18]). The solution of Eq. 01 
can be written in a form similar to Eq. QI8p . namely 



M{p) 



exp(-p) 

P 2 Pr{p) 



/3 r (u)R f {u)du. (43) 



The main feature of this solution is the divergence of the 
integral in the right-hand side as p — > 0, owing to the 
square-law dependence of the variance on x at large x. 
This singularity should be treated with care. 

To do this we note that F(p) oc a.f$p~ 2 at small p 
and hence the difference 2p 2 F 2 (p) — 2a 2 p~ 2 is regular 
at p — > 0. Therefore, it is convenient to define an entire 
function n{p) = (3(p)Rf(p) — 2a 2 p~ 2 . In terms of this 
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function the solution M(p) can be expressed as follows: 
cxp(-p) 



M(p) 



P 2 P(p) 



2a i° r ( v. 

— 1" fc/,o — / K(u)du 

P Jo 



(44) 



The asymptotic expansion of M(p) near its pole of the 
third order is of the form 

M( P ) = ^ + kf -" + ^ + (45) 
pi p 

Applying the inverse Laplace transformation, we find the 
asymptotic form of v(x): 

v(x) = a 2 x 2 + (kffi + 2a 2 j )x + fc/ j0 + «/. + 1 ; (46) 

with an exponentially small error term. Subtracting 
f 2 (x), we find v(x) — f 2 (x) = fj, r (x + 1) where fi r = 
kfo + 2a fo is the specific variance of the fill factor (at 
x — > oo), 



Mr = 3a />0 + 2 / /3 r (u)p^(u)(e p(z - 1) + 
Jo 



+2 



d«. (47) 



Integrating by parts the second term and rearranging the 
result, we rewrite \i r in the form 



fl r = 2 



a f (u) 



(e~ u + e~ ul - 2e- u(l+1 ^ du 



-2 



a2 M e - U (i+i) 

o (3 r {u)u 2 



{(l + l)u 



e~ u + e~ ul 



where 



a f (u) = a/,o 



Pr(y)dy. 



- 2) du 

(48) 

(49) 



In the limit 1 = 1, Eq. (|48|) reduces to the standard 
RSA result first obtained for a lattice RSA model by 
Mackenzie [2l| . The numerical value of the Mackenzie 
constant //o = 0.0381564.. corresponds to the Fano fac- 
tor $ = 0.0510387.. (see Q for detailed estimates of the 
variance). Expression (|4"8")) for shrinking particles has the 
same structure as the corresponding formula for the stan- 
dard RSA model (fixed-size CPP). 

Specific variance \i r of the filling factor obtained from 
Eq. (|48|) is plotted in Fig. 7 against the excess parking 
length factor I. We also plot the ratio of the filling factor 
variance to its mean value u r ja,f.Q, known in the high- 
energy detector physics [IcJIlIl as the Fano factor $. We 
see that while the fill factor variance decreases with I (as 
does the mean fill factor itself, cf. Fig. 3), their ratio $ 
increases with I due to the steeper decrease of the mean 
fill factor. 
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FIG. 7: Specific variance and the Fano factor (shown by 
curves with primes) as functions of the shrinking ratio I. 
Shown are results of exact calculations: (1,1') - scaled re- 
sults for particles shrinking to center, (2,2')- results for parti- 
cles randomly shrinking to the ends [based on Eq. (|48[) ]. and 
(3,3')- results of solution neglecting edge effects; also show re- 
sults of approximation for f(x) and u(x) — f(x) 2 at x > I + 1 
allowing edge effect (4,4') and without it (5,5'). 



Similar analysis of the model in which particles shrink 
to the center shows that for this case the variance asymp- 
totic is given by a constant [i c = 2/(1 + l)/xo, so that it 
decreases as fast as fill factor. As a result, the Fano fac- 
tor is not modified in this model. This is in line with 
the fact that both the number of adsorbed particles and 
their distribution is identical to that for a standard RSA 
with the modified particle size. 

Figure 7 also shows the variance and the Fano factor 
obtained by an approximate solution of Eq. @ based 
on linear approximations for both the mean waste length 

ix + d for 



w(x) 



- b and the variance u(x) — w(xY 



x > 1 + I recursively combined with exact expressions in 
the region x < l + l. We see that the difference between 
the approximate and exact solutions is quite substan- 
tial. Moreover, the approximate solution exaggerates the 
growth of the Fano factor with I. 

Expanding particles. In this case it is convenient 
to solve equation pip for the gap (waste) variance, us- 
ing Eq. (f2"8"|) for the boundary conditions. We define a 
function 



N(p) 



- pX 



u(x)dx 



(50) 



that satisfies an equation of the same type as Eq. (|32[) . 
in which one should replace G w (p) — > R w (p), where 



R w {p) = 2^M 6(op) + 2 p 2 (f w ( P ) + ^Map) 

pZ I pZ 



with 



J 2 (ap) 



re^dt. 



(51) 
(52) 
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1.4 1.6 U 

1 (enlargement ratio) 

FIG. 8: Specific variance (dashed curves) and the Fano factor 
(solid curves with primed labels) as functions of the expan- 
sion ratio I = 1/a. Exact calculations with Eqs. (48-50) are 
shown by curves (1,1'), curves (2,2') correspond to edge effect 
neglected, and curves (3,3') represent the linear approxima- 
tion of w(x) and u{x) — w(x) 2 at x > I + 1. 



Similarly to the standard RSA case, N(j>) has a third- 
order pole, i.e. at p — > 

.TV ^ a ^-0 , (k w ,e + a w Q ) 2 k we +a wQ 

N (p) = + — ' — 2 — — H — • ( 53 ) 



P 



P 



Separating out the regular part of N(p) and rearrang- 
ing the terms, one obtains the constant k w>e in the form 
k w e = K\ + K 2 + K 3 , where 



A'i = 2 



Ko = 4 



0(u) 



(b{au)J 2 (au) + J 2 (au)) du, (54) 



1 - 2e u b(au)) Ji(au)du, 



(55) 



and 



= -4 / " „ e~" ((e"" - l)6(o«) + u) du. 



o (3(u)u 2 



Here 



a w (u) = a w _ - 2 



o P 



(56) 



b(ap)J\{ap)dp. (57) 



Variance of the gap distribution is obtained by the in- 
verse Laplace transformation of N(p). At large x one has 
v(x) — w 2 (x) = k Wt e(x+ 1). The constant k w ,e describes, 
with exponentially small errors, the linear dependencies 
of both the filled length and the gap variance. 

With the substitution a w (u) = (u + l){3(u) — af(u) and 
integration by parts, formulae <|54ti56[) can be rearranged 
after some lengthy algebra into a form similar in structure 
to Eq. P5|) . However, additional terms at a ^ 1 make it 



rather unwieldy. Both forms arc equivalent for numerical 
integration. 

Variance of the gap distribution calculated with Eqs. 
(|54j|56[) equals the filled length variance. Figure 8 dis- 
plays specific variance along with the Fano factor as func- 
tions of I = 1/a. Note the increase of both the variance 
and the Fano factor at small I. This increase is appar- 
ently due to the singular contribution to OPDF from edge 
effect. For comparison, we show the v(x) and $ calcu- 
lated neglecting the edge effect (this is accomplished by 
taking b(ap) = 1): no increase is seen in this approxi- 
mation. Both the variance and the Fano factor vanish 
with increasing I owing to the fact that variably growing 
particles provide tighter filling. 

Figure 8 also shows the variance and the Fano factor 
obtained by an approximate solution of Eq. (|11[) based on 
the linear approximations for both the mean waste length 
w(x) — ax + b and the variance u(x) — w 2 {x) = cx + d for 
x > 1 + 1 combined with exact expressions for x < 1 + I . 
In contrast to the case of shrinking particles, the differ- 
ence between the approximate and the exact solutions 
becomes smaller at larger expansion ratios, since both 
solutions predict a rapid decrease of the Fano factor with 
I. Note that in this case the overall variation is well de- 
scribed by the approximate solution. 



V. CONCLUSIONS 

We have considered a generalized 1-dimcnsional ran- 
dom sequential adsorption problem, where particles 
shrink or expand upon adsorption. Using a recursive ap- 
proach, we obtained exact analytical expressions for the 
filling factor and its variance. 

In the model where particles shrink to their centers, 
we find that both the filling factor and its variance scale 
as the recursion period, so that the Fano factor remains 
unchanged. But in another model of shrinking, where 
particles shrink by randomly retracting their endpoints, 
the Fano factor increases with larger shrinkage ratio, due 
to the weaker decrease of the variance compared to that 
of the filling factor. In the case of adsorption of expand- 
ing particles (within the allowed intervals), we find that 
the variance first increases with the expansion ratio I and 
then decreases, showing a maximum at relatively small 
values of I. These nontrivial results are due to edge ef- 
fects, resulting in a modified one-particle gap distribution 
function, p(y), near the edges of the adsorbing interval. 

The developed approach can be applied to other mod- 
els, so long as they correspond to a piccewise constant 
p(y). The results are shown to have exponentially small 
corrections for finite adsorption length x; for x > 7 I they 
become practically exact. 

The analytic results have been compared to approx- 
imate solutions based on a method frequently used in 
semiconductor physics. The approximate solution is rea- 
sonably accurate in estimating the fill factor for cither 
shrinking or expanding particles. However, the same 
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method is less accurate in estimates of the fill factor vari- 
ance. In the model where particles shrink by retracting 
their endpoints the approximate method overestimates 
the effects of shrinkage on both the variance and the Fano 
factor. We also assessed another common approximate 
method, based on the mean final energy distribution ap- 
proach, and found that it gives only qualitative trends. 

Our quantitative results have important applications 
to the problem of energy branching in high-energy parti- 
cle propagation through a semiconductor crystal, where 
the model of shrinking particles (where particles shrink 



symmetrically about their centers) naturally accounts 
for the fact that the impact ionization threshold is 
larger than the energy gap for electron-hole pair gen- 
eration. The alternative model of shrinking particles, 
where shrinking occurs by retracting one of the particle 
endpoints, further accounts for the decreasing density of 
states at low particle energies. 
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APPENDIX A: RESTRICTED EXPANSION 
MODEL, COMPARISON WITH RATE 
EQUATION RESULTS 

In the model of RSA with expanding particles stud- 
ied in H, the authors considered the case of restricted 
growth, when the deposited particle begins to grow only 
if the interval in which it is adsorbed is larger than its 
final size. For this model, Eq. (jlip remains valid for in- 
tervals x > 1, but the initial conditions to Eq. (TTTj) for 
w(x) (at x < 1) are of the form 

!x, < x < a 

x — a, a < x < 1 (Al) 
x-1, 1 < x < 1 + a 

To solve this model, we shall use the same technique as 
in Section IIII1 so as to compare the results of the kinetic 
and the recursive approaches. To do this, we multiply 
Eq. (|11[) by x — a, make the substitution x — > x + 1 and 
then take Laplace transformation. We obtain an equation 
similar in form to Eq. (|30[) where, however, J(p) must 
be calculated using the revised initial conditions of Eq. 
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(TAT]) . This results in 



J{p) = / w(x)e- p: 
Jo 



1 



= -[j 1 (ap) + e~ ap J 1 {{l-a)p)]. 

(A2) 

The solution of Eq. (|30j) and hence a^o are of the form 
similar to Eqs. ([55)) and ([55]) in which one must use 
G w {p) defined by Eq. (|3~4"| but with J\{ap) replaced by 
the sum of two terms in the square brackets of Eq. (|A2[) . 
Calculating the filling factor as ctffi — 1 — &w,Qi we nn d 



a/.o = 2 / I3{p)b(ap) [ae 
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+ (1 - a)e~ p ] dp. (A3) 



In Ref. Q, the contribution of expanded and non- 
expanded particles was calculated separately for a finite 
time t, using the rate equation approach. In order to 
make the comparison, one has to take an appropriate 
linear combination (Eq. (15) and Eqs. (31), (32) of Ref. 
0) in the limit of infinite time. This yields (in the nota- 
tion of 0]) the fill factor 8 in the form: 



8 = f°° F(t) {l + {a - 1) [1 + (a - l)t] e"^ 1 )*} 



dt. 

(A4) 

To compare Eqs. (|A3|) and Eq. (|A4[) , we multiply the 
right-hand side integrand in (|A3[) by unity in the form 
1 — e p + e p . Then, the contribution /, proportional to 



1 — e p , can be rewritten as 

I = J d (/3(;p)) p \ae (1 - a)p + (1 - a)] (A5) 

and then integrated by parts to give 

/ = - J p{p) {a [1 + (1 - a)p\ e (1 ~ a)p + 1 - a} dp. 

(A6) 

The remaining part iY is given by 

f°° 

11= (3{p) ae (1 - Q)p + (1 - a) [2 + (1 - a)p] dp. 
Jo L - 1 

(A7) 

The snm I -\- II gives af^o in the form 

poo 

a f , = / Hp) ae (1_a)p + (1 - <0(1 + (1 - a)p dp. 
Jo L J 

(A8) 

Equation (|A8|) . after the substitutions ${p) = 
F(p) exp[— (1 — a)p], a = 1/a and p = at, transforms 
into an expression identical to the 8 given by Eq. (|A4[) 
derived from 

It should be noted, however, that while the fill factors 
are identical, success of the kinetic approach is hard to 
extend to variance calculations, where it appears to be 
much less effective (see e.g. Q). 



